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Abstract 

We estimate the properties of clouds that form as a result of the inter- 
action of the surface of a disc embedded in a supersonic outflow and heated 
by external irradiation. We flnd two regions where short-lived clouds are 
injected into the outflow and accelerated by ram pressure. We inden- 
tify these with a broad line region and an intermediate line region. We 
compute the line strengths and profiles from a model of the cloud dynam- 
ics and compare the results with the observations of the nearest Seyfert 
galaxy NGC 4151, which, in its various luminosity states, has provided a 
large amount of emission line data. We show that the model parameters 
for this galaxy are constrained by just the CIV A1549 line profile in the 
high luminosity state of the nucleus. We review briefly some of the data 
on broad line ratios, proflles, profile variations and transfer functions and 
show how much (if not all) of this data can be accounted for in the model. 
We also indicate how the X-ray data can be fitted into this picture. 

1 Introduction 

We can distinguish various approaches to understanding the broad emission lines 
from active galactic nuclei. In phenomenological models we attempt to obtain 
a distribution of photoionised gas that will reproduce the main features of the 
line ratios and profiles without regard to the dynamics of the gas or its relation 
to the system as a whole. Early investigations resulted in the now universally 
accepted cloud picture, but were ambiguous with regard to the cloud kinematics. 
More recent investigations of time lags between continuum and line variations 
restrict the scale of the cloud distribution and point to complex kinematics and 
a stratification in cloud properties with radial distance from the ionising source. 
To accomodate these observations, phenomenological models must incorporate 
a range of parameters or arbitrary functions (for example, Robinson 1995). An 
alternative approach is to try to model the cloud system dynamically. Gravity, 
magnetic fields, ram pressure and radiation pressure have all been considered 
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for the acceleration of pre-existing clouds (Carroll & Kwan 1985, Emmering et 
al. 1991, Capriotti, Foltz & Byard 1979, Shields 1978, Blumenthal & Mathews 
1975, Mathews 1993). Two phase equilibria and shocks have been considered for 
the formation process (Mathews & Doane 1990, Perry & Dyson 1985). Models 
of this sort have been used to address either individual galaxies or classes of 
galaxies (CoUin-Souffrin et al. 1988, Terlevich 1994). 

The present paper stems from a phenomenological model (Cassidy & Raine, 
1993) in which clouds are supposed to be injected from the surface of an accre- 
tion disc into a supercritical, mildly relativistic, outflowing spherical wind where 
they are accelerated and destroyed. Although the existence of such a wind is an 
assumption for the formation of which there is no firm theoretical evidence, we 
showed that this picture gives a good account of the large variety of line pro- 
file shapes that occur in active nuclei and we were able to relate characteristic 
features of the profile types to aspects of the physical model. This work was in 
turn based on a theoretical model (Smith & Raine 1984) of the formation and 
acceleration processes in a thermally driven wind from the disc surface which 
we now know to be untenable for the broad lines (Sect. 2). Here we present a 
new physical picture of the interaction of an accretion disc in an ambient wind 
in which the broad line clouds are formed as a result of Kelvin-Helmholtz insta- 
bilities at the interface. A further new feature of the model is that the thermally 
driven outflow from the disc provides a region of line emitting clouds (the inter- 
mediate line region, ILR) with properties intermediate between the true broad 
line region (BLR) and the narrow line region (NLR). The theory supports the 
phenomenology described in Cassidy and Raine. It enables us to go some way 
towards a detailed description of the broad line emission from NGC 4151, one 
of the most enigmatic (and hence theoretically challenging) of galaxies. In Fig.l 
we give a schematic diagram of the main features of the model. 

The picture we present attempts to address seriously from a theoretical view- 
point the wealth of detailed observational data that point to a complex, strat- 
ified, emission line region. Our model is similar to Emmering et al. (1991) 
except that they take magnetic forces to be partly responsible for the injec- 
tion and acceleration of cloud material from the surface of the disc into the 
broad line region. It differs from biconical models (Zheng, Binette & Sulentic 
1990) in that our cloud outflow fills a conical shell rather than the interior of 
a (double) cone and our cloud motion inherits an azimuthal component from 
the disc. Conversely, we differ from pure disc models (e.g. Shields 1978) since 
in our case there is a radial component to the cloud velocity. In contrast to 
the two component model of Collin-Souffrin et al (1988), we obtain also the low 
ionisation lines from the cloud component and not from disc emission. The two 
component model is certainly not applicable to NGC 4151 in any obvious way 
since in this source the CIV A1549 and Mgll A2798 line profiles are similar to 
each other and both different from CIII] A1909. 

Clearly our picture is not spherically symmetric so differs from, for example, 
Dyson and Perry (1985), Mathews (1993), Terlevich et al (1994) and many 
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others. Note that we do inject short-hved clouds into an outflowing wind but 
(in contrast to, for example, Dyson and Perry) the cloud material does not 
originate in the wind. Nevertheless, wc agree with many authors in taking the 
components of the active nucleus to be a disc, a nuclear source of both radiation 
and particle flux and stars. We difi'er in that we explore the interaction of the 
nuclear wind and the disc as the source of emission line clouds. Of course, 
our nuclear wind must interact with the stars in the BLR as in the Dyson and 
Perry picture; but we assume that, for NGC 4151 at least, cloud formation from 
shocks induced in the wind by the present population of supernova remnants is 
too small to make a signiflcant contribution to the BLR (although it may make 
some contribution to the ILR). We also neglect the contribution from bloated 
stars (Alexander & Netzer, 1994) and the stellar ionising continuum. 

NGC 4151 exhibits changes of state that can be described as a sequence 
of intermediate Seyfert types. Ulrich et al. (1991) distinguish between a high 
luminosity state (fAi455 > 15 x 10""'^'' erg cm^^ s^^ ^^^) SiJid a low luminosity 
state (fAi455 < 3 X 10"^'' erg cm~^ s~^ A~^). At its highest luminosity NGC 
4151 shows features characteristic of Seyfert I nuclei with broad optical and UV 
hnes (up to 30 000 km s"! FWZI) and a high H/3/[0III] ratio (Clavel et al. 
1990). At its lowest luminosity the broad lines have all but disappeared and 
the nucleus is described as a Seyfert 1.9. Changes between states take place 
on a time scale of weeks to months, short even in comparison with other such 
intermediate Seyferts. However, the extreme luminosities occur very rarely. In 
section 3 we shall show that the CIV line profile in the high luminosity state fixes 
the parameters of the model. We consider also the likely effect of reductions of 
luminosity using the theory of section 2 as a guide. In section 4 we explore the 
extent to which the model might yield agreement with the extensive range of 
observational data on this galaxy. 

2 The Theoretical Basis 
2.1 Introduction 

In this section wc shall show that for an accretion disc embedded in a supersonic 
outfiow, and illuminated by an external radiation field, there are two regions 
from which clouds are injected from the disc surface into the ambient outflow. 
The inner of these regions yields clouds with the properties of the broad line gas 
(BLR) and the outer region can be identified with the intermediate emission 
line region (ILR). We shall anticipate this identification and refer to the two 
zones of the disc as the BLR and ILR. The theory will be developed in general 
and applied to the source NGC 4151; thus we shall write some numerical factors 
in equations in terms of the parameters for this galaxy, even though, of course, 
this will require us to assume the results of the fit to the observations in sections 
3. 
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External illumination by X-rays raises the temperature of material at the 
surface of the disc to the Compton temperature. At sufficiently large radii 
this material is then hot enough to escape from the gravitational pull of the 
black hole (Begelman & McKee 1983). We have shown that if this material is 
subject to the ram pressure from a nuclear wind then unstable cooling leads to 
cloud formation (Smith & Raine 1984). Our original proposal that this is the 
mechanism for the formation of broad line clouds is, however, untenable. For a 
typical Compton temperature of the incident radiation field, — 10^, mass loss 
from the disc can occur only at radii greater than Tc — lO^^mg/Ty cm for a black 
hole of mass lO^mgMQ, where T7 — lO^Tc Even if the Compton temperature 
of the radiation incident on the disc differs from that in the line of sight to the 
observer (which is not the case in our model of NGC 4151) the inferred radius of 
the supposed BLR would still be larger than permitted by variability arguments. 
Furthermore, if the Compton temperature is lO^K, the density in the clouds, at 
temperature lO'^K, will be of order 10* cm""^ (see equation (||)), too low to be 
compatible with the typical inferred value of the BLR ionisation parameter of 
logU ~ — 2. In addition, in order to obtain logarithmic line profiles in the line 
wings, we assumed that the clouds are long-lived. The longevity of the clouds is 
inconsistent with acceleration by a supersonic wind (since the clouds loose mass 
through the flanks on the order of a sound crossing time). Short-lived clouds in 
this model cannot produce acceptable line profiles (Cassidy, 1994) nor provide 
the variety of non-logarithmic profiles now observed in other sources (Robinson 
1995). The requirements are relaxed somewhat if, as may be the case in NGC 
4151, the Compton temperature is as high as lO^K, but the overall conclusions 
are unaffected. 

Nevertheless, this mass loss must occur, and, if the nuclear wind is present, 
cloud formation is a consequence. We shall show that these clouds must con- 
tribute to an intermediate line region. In the low luminosity states of NGC 
4151 the ILR dominates the line emission; from this we shall deduce that these 
clouds may be sufficient to provide the main component of the ILR emission in 
this source. 

We shall show further that within the Compton radius, rc, the illuminated 
surface of the disc in the presence of the nuclear outflow is Kelvin-Hclmholtz 
unstable. These instabilities, we argue, will grow and be entrained in the wind 
as BLR clouds. 

In Cassidy & Raine (1993) we used a phenomenological picture based on the 
injection of clouds from the surface of the disc by an unspecified mechanism to 
generate the broad line region. This work gave a range of line profiles that could 
be matched to the variety of profile classes observed, but at the expense of a 
large number of adjustable parameters since the physical basis of the model was 
not specified. In the present paper we take a different approach: we construct 
a theory of the cloud injection process which enables us to relate some of the 
parameters of the system by self-consistency arguments. We then apply the 
results to one of the most difficult examples, NGC 4151, for which we can find 
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a set of values for the remaining parameters from a single observation of the 
CIV A1549 line profile at its peak luminosity. In particular we fix the angle of 
the disc to the line of sight. In a later paper wc shall show that NGC 4151 (in 
each of its various luminosity states) can indeed be considered a typical active 
nucleus viewed from the particular line of sight, just over the top of the disc, at 
which the nearside of the flattened distribution of BLR clouds appears directly 
in front of the central continuum source. 

Thus, to summarise, the active nucleus is supposed to contain a black hole, 
of mass TOg in units of IO^Mq; there is a disc, which is taken to be a source 
of material, but about which we make no assumptions other than that it is 
opaque with non-zero albedo; we require a central source of luminosity Lj., 
parametcriscd by Ir = Lr/Ls, where Lb = 1.3 x W^^ms is the Eddington limit 
luminosity; and we have a mass outflow which carries a flux of kinetic energy 
Lm- This matter flux, the Eddington luminosity and the disc inclination are 
the only true free parameters in the theory, apart, that is, from dimensionless 
parameters of order unity, behind which we hide our ignorance of some of the 
details of the physics. The results depend significantly on only one of these 
dimensionless parameters, the disc albedo. We also have to make assumptions 
about the continuum spectrum in those bands where it is not measured, but 
our results are not sensitive to this choice, except for some dependence on the 
Compton temperature of the disc illumination (which is not necessarily the 
same as the observed Compton temperature). As noted below, we prefer a 
Compton temperature of the illuminating radiation closer to 10* K than 10^ K. 
The wind must be sufficiently optically thin to be compatible with the X-ray 
variability. We use the observations of NGC 4151 to calibrate this constraint. 
The interactions in this system then lead to a scattering of radiation in the wind 
on to the disc, a zone in which the disc surface is Kclvin-Helmholtz unstable, 
which is the BLR, and a zone in which the disc surface is the source of a 
thermally driven wind, which forms the ILR (Fig.l). 

2.2 Wind Properties 

The supersonic nuclear wind is taken to be spherical with a constant outflow 
velocity Vyj = 10^°i;io cm s~^ and to carry a power L„ in bulk motion. The 
wind density, n^, is then given by 

= 4 X 10\Lm/LE){m8/0.3)r^eV^o^ cm-^ (1) 

and the ram pressure is 

P,n = 7 X 10\L^/LE)ims/0.3)r^e^v^o^ erg cm'^. (2) 

The wind is fully ionized by the radiation field and heated to the Compton 
temperature out to a radius rad = 3 x cm. Beyond this adiabatic 

cooling will start to become important but the wind temperature will not drop 
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substantially until r = IQ^rad (Smith and Raine 1985). The Mach number of 
the outflow is of order lOOwioTg"^^^ and the wind is typically hypersonic. The 
optical depth to electron scattering is 

= 2 X 10-2(m8/0.3)(i„/LE)rr6'vro^. (3) 

So, for NGC 4151, provided Lm/LE < 5, we have Tu, ~ 1 at r = 10^^ cm which is 

consistent with large scale X-ray variability on times of order 10"'' s. (Yaqoob ct 
al. 1991, 1992). Then, consistently, = 10^" cm s~^ » v^sc at 10^^ cm (since 
Vesc = 10^°((?n8/0.3)/ri4) 2 cm s~^). We cannot reduce the outflow velocity 
significantly below t'lo = 1 for a high value of Lm/ Le, as will be required to fit 
the high luminosity line profiles below. The nuclear outflow must therefore be 
mildly relativistic for consistency. 

2.3 Properties of the disc 

The radiation falling on the disc will significantly affect the surface if the energy 
falling on the disc is greater than that produced locally by viscosity in the disc. 
In particular the surface layers will be heated to the Compton temperature. Note 
that under these circumstances the disc may still have an optically thick, cool, 
geometrically thin core (Mobasher & Raine 1990) at this radius. Suppose the 
disc is illuminated directly (rather than by scattered radiation). Let the surface 
be given by z = z[r). Then the luminosity intercepted by the disc between radii 
r and r + dr is proportional to 2'iTrdz/r'^. Let the flux falling on the disc be 
{(pLr/4:Trr'^)dz/dr between radii ri (where = 1) and Vmax- (This defines (j).) 
Taking z{r) to be given by a scale height in the hot surface layers of the disc 
when external illumination dominates gives z{r) = xIO^^T^^^Tiq^ {ms/0.3)~^^'^ 
cm. For the local dissipation in the disc it is reasonable to assume a steady 
state. But the irradiation from the central source may change on short time 
scales. In terms of a time-averaged luminosity parameter 1^, we therefore have 
Lr = 0.1{lr/lr)Mc^ for an accretion eSiciency of 0.1 and an accretion rate M. 
This gives the radius at which external illumination equals the local rate of 
viscous energy generation as 

nn ~ 10^^(^)^/^(^)^/^T3-V3(|l)V3 enr. (4) 

This therefore gives an estimate of the inner radius of the BLR which conse- 
quently ranges between about 10^^ cm and 5 x 10^*' cm as lr changes by a 
factor 10. The outer radius of the BLR depends on the flaring of the disc in the 
presence of the wind and is calculated in section (2.7) below. 

The surface temperature of the disc is the Compton temperature. Beyond 
a radius rc = 3 x 10^^(TO8/0.3)Tg"^ cm this exceeds the temperature at which 
material is gravitationally bound (Begelman et al 1983). The material escaping 
from the disc has been raised to a pressure P = Pj-adf^'c min where S is the 
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ionisation parameter S = Prad / Pgas and is its value at which the gas 

can just remain in the hot phase. The escaping material therefore has a density 
jO = P/Cj . Exposed to the ram pressure of the nuclear wind this material cools 
under constant pressure to lO^K. Using = 3 as a typical value for a 

Seyfert spectrum, the density of the clouds becomes 

These clouds constitute the ILR. Begelman et al. (1983) find that there is a 

— —1/2 

fully developped disc wind extending from rc to a radius Viso — 3(Zr/0.1)Tg Tc 
within which material is rapidly heated to Tc. Thus the ILR extends from within 
3 X 10^'' cm to beyond 10^^ cm. 

Most models of accretion discs are optically thick to the continuum in the 
BLR, but not usually in the ILR (CoUin-Souffrin 1987, Mobasher & Raine 1990). 
Therefore the disc obscures any emitting material in the BLR on the far side 
from the observer, but not in the ILR. In the BLR we expect the albedo of the 
disc to be significant at least for resonance lines. Thus, we shall have to take 
account of reflection from the disc surface. 



2.4 Disc-wind interaction in the BLR 

Consider now the interaction between the disc and the nuclear wind in the BLR 
region. Suppose, contrary to our discussion above, that the disc surface were 
to remain cold and therefore essentially unexposed to the ram pressure of the 
wind. Then the pressure in the disc is balanced by the thermal pressure, n^kTc 
in the wind. If a fraction (j) of the radiation field falls on the disc, the ionization 
parameter (Begelman, McKee & Shields 1983) is therefore 

S = ^ = /^^2xl0^^ (6) 

for an inverse Compton temperature Tc = lO^Tg K. By running the CLOUDY 
photoionisation code for an average quasar spectrum (Mathews and Ferland 
1987), a power law with an energy index of -0.9, and a power law with the 
canonical slope of -0.7, we find that the material must be hot, with temperatures 
of 1.5 X lO"", 1.0 X 10^ and 1.6 x 10^ at values of S greater than of order 10^, 10 
and 1 respectively. The value of S from equation (^) is therefore much greater 
than the maximum for the cold phase. Thus, contrary to our assumption 

(and in agreement with section (2.3)) the disc surface cannot remain cold. But, 
also, the disc surface cannot be hot and be exposed to the full ram pressure of 
the wind since in that case 

S=^ = (|^)(^), (7) 
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smaller than a critical value 2c = S'^ the minimum for the hot phase. So, 
we can assume that the disc surface adjusts itself so that the disc is exposed to 
a pressure fPm just insufficient to maintain the cold phase: 

f = c^lr{^)C-T)^- (8) 

If, for example, Lj^/Le < 5 and lr(j) = 0.01, we have / > 2 x 10~''(3/Sc). 

From this we can deduce the density in the surface of the disc, n^, since the 
disc pressure must balance /P,„ at the Compton temperature. We obtain 

= l.6xlO^{(t)/0.1){lr/0.1){ms/0.3)E^\^eT^^cm'^ (9) 

Next we show that this arrangement is Kclvin-Hclmholtz unstable at the 
disc-wind interface. Rotation and infall of the disc material can be neglected 
since the rotational and infall speeds are small compared with the wind speed. 
Thus, for Ud ^ Uw the instability operates on wavenumbers (Chandrasekhar 
1961) 

k > (10) 

where g is the gravity at the surface. Suppose the surface is at height z = , 
where H is the pressure scale height at the temperature of the disc atmosphere, 
lO^Tg K. Then, the vertical component of the acceleration due to gravity is, 
approximately for z ^ r, 

9 = - 6.4(™8/0.3)T|rr/C cms'^. (11) 

The height C is determined by pressure balance in the disc. For any disc 
model it is unlikely that C ^ 1 (because the pressure falls off rapidly with 
height as exp(-C^)), or <C 1 (because the external pressure does not dominate 
the internal disc pressure). It is therefore sufficient to take C = 1- 

Using equation (|^) for the density in the disc, we obtain 

k > krmn = 2.6xlO-i6(m8/O.3)(i£/i„,)(0/O-l)(/r/O.l)Cfior8''/'r-//'s-icm-^ 

(12) 

and potentially large columns 

uXTh = 4 X 10^^{L„,/LE){ms/0.3)r-^^^\^o'Ts'^\'^ cm-^, (13) 

where 

XTh = 27r/fc„„„ = 2 X lOi«(O.Ol/0Z,)(O.3/m8)(im/iB)Scfro'rV-^/'T8 cm 

(14) 
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is the length scale of the growing perturbations. Note that at the inner radius 
this is larger than the BLR, so this maximum, which is calculated assuming 
a uniform system, cannot be attained. The outer radius of the BLR will turn 
out to be of order 10^^ cm, (section 2.7) so ^KWi^out) ~ ^out i.e. the maximum 
scale of the growing mode is of the order of the size of the BLR. However, 
the growth time decreases for smaller wavelengths (and is infinite for the largest 
mode (equation (17)) so the disc surface will develop large rising rings of hot gas 
on various scales within the BLR. This gas is exposed to an increased pressure 
from the impact of the nuclear wind at which the material cannot remain hot. 

The main contribution to cooling above about lO^K is thermal bremsstrahlung, 
so the cooling timescale in the disc surface is 

to ~ 10-27^1^1/2 10'(0.1/</.)(0.1/«.)(0.3/m8)r|/2Ser?6 sec. (15) 

This is to be compared with the sound-crossing time, tg. For material of column 
density, Ncoi, we have tg = Ncoi/nCg and hence 

ts/tc ~ {Ncoi/lQl'^)T^\ (16) 

independent of any parameters of the system. Thus, the hot bubble can cool 
coherently only in a relatively small column of density about 10^^ cm~^. Cooling 
to 10^ K will tend to be one- dimensional producing a surface layer of this column 
density. The growth time of the instability is 



tg = { 7 "' )-^(^-l)-i (17) 



which, for scales Xkh <S ^xif depcnids only on column density so is unaffected 
by cooling. Perturbations on the cold surface grow on a timescale of order 
10~Hs. Thus, the surface will tend to ripple as it cools into clouds of number 
density, assuming cooling at constant pressure, 

3 X 10^2 , 6 ^, Ir ^, 1 ^ o 

rfg 0.1 0.1 Tg^c 

Where the timescale for cloud formation, tc, exceeds the variability timescale 
the factor of Ir in this equation should be replaced by its time averaged value 
Ir- But for most, if not all, of the BLR we can ignore changes in cloud density. 
These clouds will be entrained into the nuclear wind. 

We can now summarise how the instabilities develop. Hot rings grow on the 
disc surface and start to rise into the wind. The material at the surface of the 
rising bubble is exposed to a higher pressure and a shell of column of Nc — 10^^ 
cm~^ starts to cool towards lO^K. The cool material is still Kelvin-Helmholtz 
unstable on scales less than Nc so cooling does not suppress the growth but will 
tend to break the rings up into clouds.. Thus, the effect of the instability is to 
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inject a spray of clouds of column density 10^'^ cm^^ into the nuclear wind. Let 
Ac = Nc/rid, the thickness of the shell of disc material that cools in a sound 
crossing time. The efficiency with which the hot bubble is injected as cold clouds 
depends on the ratio Xc/Xkh, since this is the fraction of material in the bubble 
that cools to form the surface shell. Thus the number of clouds injected as a 
function of radius is: 

^c/XTh = 3 X m-\LE/L^){0.3/ms)vioT^/\rli\ (19) 

The important point is that this increases with radius as r^/^ thereby weighting 
the outer clouds. 

The density of the cold material at lO^K is a factor lO^Tg greater than the 
surface density of the hot disc (equation (|l8|)). Such densities (of order 3 x 10^^ 
cm~^ for the innermost clouds in NGC 4151, falling to 3 x 10^° cm~'^ in the 
outer BLR) are much higher than usually attributed to the BLR. We shall find 
such densities to be precisely those required for agreement with the CIV/CIII] 
line ratios and profiles. 

2.5 Cloud injection 

We argue next that these cool regions are injected into the nuclear wind with 
some velocity Vi and accelerated to some velocity Vc- Consider an elevated 
region of surface wave. The growing protuberance provides an obstacle around 
which the wind must flow. While the face of the perturbation is exposed to 
some fraction of the ram pressure, the top surface will be subject to a smaller 
pressure closer to the thermal pressure, Pw, in the wind. Since the disc material 
is on average at fP,n ^ Pw the incipient hot bubble, area A, is accelerated 
upwards to a velocity Vi given by 

fP„A = {A\KHndmp){^). (20) 

This gives vf — Pd/rapUd or Vi = Cd, the sound speed in the hot disc surface. 
The sprays of cold bubbles are injected with velocities normal to the disc of 
order lO^Tg cm s^^. Note that, even if the Compton temperature of the disc 
surface is of order of the observed average of lO^K, mechanical heating of the 
atmosphere may raise the temperature above Tg ~ .1 so injection velocities 
above 3 x 10^ cm s~^ may be possible (King & Czerny 1989). (In fact we 
shall take Tg ~ 1 for NGC 4151 on other grounds; see sect. 3) The clouds will 
also inherit a rotational velocity from the disc. The disc surface must flare to 
a scale height at the Compton temperature at the outer edge of the BLR (ie 
H/r ~ 0.5) so the rotational velocity will not be Keplerian. However, since 
H/r < 1 throughout the BLR and H/r < 0.1 for much of the region, we shall 
take the azimuthal component of the cloud velocity to be Keplerian on injection. 
In fact, the drag on the cloud as it moves across the radially flowing wind reduces 
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the rotational component v,/, significantly in a time (cs/w</))ts, where ts is here 
the sound crossing time in the cloud. Since the azimuthal velocity in the disc 
is supersonic, the rotational component is strongly attenuated. So although it 
cannot be neglected in forming the line profiles (since it is largest in the initial 
part of their motion where the clouds spend the most time) departures from 
Kepler velocities are unlikely to be of crucial significance. We therefore assume 
the clouds inherit a Keplerian velocity from the disc on injection. The injected 
clouds will be subject to the full ram pressure of the wind. This will act to 
compress and accelerate the clouds but also to destroy them. We expect the 
cloud survival time to be of order a sound crossing time in the uncompressed 
cloud. Thus, approximately, for the cloud acceleration, wc have (ignoring mass 
loss) 

PcAXc^^p^vlA (21) 
or, taking t = ts = Xc/cg (c^ is the sound speed in the cold cloud), 

Ve^CsC-f){-)-r'cs, (22) 



Pc 



or 



{(j)/.l)[lr/.l)Vl0 Le 

using equation (8). For / ~ 5 x 10"'' we obtain speeds up to 2 x 10^ cm s~^. 
However, on the one hand, this estimate takes no account of mass loss from the 
clouds over their lifetime, whereas in the numerical computation of the cloud 
trajectories this is incorporated through a simple model of cloud disruption (see 
sect. 2.6). This estimate, of the cloud velocities should, as a result, be increased 
by a factor of, approximately, the ratio of the initial and final column densities of 
the cloud, {N{0)/N{ts)). On the other hand, the numerical computation takes 
account of gravity and the full three dimensional geometry with the overall result 
that the maximum final velocity (hence the line half-width at zero intensity) is 
of order 10^ cm s~^. 

Notice that, other things being equal, the cloud velocity (equation (p3|)) 
appears to increase as the luminosity in the radiation field, Ir, goes down. This 
is the converse of what is observed for NGC 4151. However, we do not expect the 
power in the wind to be independent of that in the radiation field. To produce 
the observed reduction in line width at low luminosity exactly in this theory we 
must require that the momentum flux in the wind, Lm/vio, be proportional to 

with (7 ~ 2. This would be the case, for example, for a pair driven wind. We 
obtain qualitative agreement provided Pm decreases with L. 

2.6 Cloud destruction 

Since the clouds are accelerated in a supersonic wind they are subject to dis- 
ruption by mass loss through the flanks, where the wind pressure is reduced 
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relative to the cloud face as the wind material accelerates round the cloud. The 
results are not too sensitive to the manner in which the clouds break up, pro- 
vided that the result of evolution is to reduce the column density, so we take a 
simple model. We choose for the mass of the cloud M{t) = M(0)(1 + t/ts)~^ 
where tg is the sound crossing time in the cloud. Thus, approximately 90% of 
the cloud is lost in two sound-crossing times. The area of the cloud facing into 
the wind and the matter density in the cloud are assumed to be constant while 
the column density decreases with the mass. 

We can give some justification for this model as follows. Mass flow through 
the flanks of a cloud at the sound speed leads, after a time t, to an increase in 
area facing the source to irr"^ = 7rro(l + Cgf/ro)^, where ro is the initial radius 
of the face of the cloud. This escaping material will be heated to the Compton 
temperature and entrained in the wind. It is therefore lost from the cloud. Thus 
the mass of the cloud is reduced by a factor {1 + t/ts)'^ while the area of cold 
cloud material remains constant. Mass must flow through the cloud to supply 
the mass loss. Assume that the density in the cloud remains constant. (This 
would be a reasonable approximation for a flattened cloud evolving in pressure 
equilibrium, which is probably not too far from being the case.) Consequently 
the column too decreases by this same factor. 



2.7 Outer limit of the BLR 

There is an outer limit to the cloud injection region set by the extent to which 
the disc surface flares to intercept the momentum flux in the wind. The disc 
pressure must be provided by new wind material further from the disc midplane 
as r increases, because the cloud acceleration at smaller r extracts momentum 
from the wind close to the disc surface. At the outer edge of the BLR clouds 
rise to a height 

Zout ~ tsVi ~ cm » H{rout)- (24) 

(0/.l)(/r-/.l) 

The momentum flux in the wind in a wedge of azimuthal extent and altitude 
= Zout /rout is, approximately 

A0 r sinOde vlcosOp^T^ ~ i(y^^^^jpll^^)l^i^^. (25) 

The maximum momentum flux in the clouds is 

A</.M„„.t;e = 2 X lQ^^logT,,M:^^^^d^Y''. (26) 

We obtain equality when ne ~ 10. Thus, for the maximum cloud injection rate, 
the outer radius of the BLR in NGC 4151 is at 10^'^ cm. 



12 



We have seen that the disc must flare to a scale height at the Compton 
temperature at the outer limit of the BLR. For NGC 4151 (mg — 0.3) the scale 

height is iJ 1.5 x VS^'^T^r% or 4.5 x 10^^ cm at Tout = lO^^cm if Tg = 1- So 
the disc flares to about .45 radians or about 30° across the BLR. Beyond the 
BLR is a region in which the hot disc atmosphere rises with z ~ if oc r"^"^ until 
this material escapes in a wind beyond 3 x 10^^ cm. 

2.8 The ILR 

The mass flow into the ILR from the thermally driven disc wind can be estimated 
by assuming the heated material flows out at the sound speed; thus 

AM = nrfCsTOp(27rrdr). (27) 

The density in the disc surface Ud is V^r'^nij^B. (equation (5))and Cs = lO^Tg^^^ 
cm s~^. Assuming, as shown in Smith and Raine (1982), that the material cools 
into clouds of column density ^ lO^'^Tg cni~^, we obtain for the sound crossing 
time in a cloud, t^, 

= 109(r/r,)2(m8/0.3)(0.1//,)(0.1/(^)T8-i s. (28) 

If the clouds survive for of the order of a sound crossing time only, we obtain 

dM = isdil = 27rlO"^Tg'/Vdr, (29) 

and hence, for the mass in clouds in the ILR, 

Mii^R = i{lrlQ.lf{m^lQ.ifT-'"^MQ (30) 

between the inner and outer radii given in section (2.3). We shall show this 
is sufficient to provide the constant narrow component of the BLR emission. 
For reference, the covering factor of the clouds (assumed spherical), is of order 
.OlTg log(3(Z^/0.1)T8"^/^) and the number of clouds is about lO^ilrclj/OmfTl'' 
These estimates all lower limits since they assume the clouds survive only a 
sound crossing time. 

3 The model profiles 

The key to the comparison of our model with observations is the fit to the high 
state profile of the CIV line. In our earlier paper (Cassidy and Raine, 1993) we 
presented such a profile fit as an example of a Df class profile. This is reproduced 
in Fig. 2 here. 

The requirement of sufficient cloud acceleration fixes Lm/ Le (equation (p3|)). 
We need Lm/LE ^ 5, an unexpectedly high value, since it implies that, in the 
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high state, most of the energy output from the central source is in the nuclear 
wind. If we require Lm/v\o oc (section 2.5) then Lm/LE ^ 1 in the low 
state and there is no overall problem in hiding the wind energy. The profile 
fit, given the column density, fixes the inner radius at rj„ = lO^^cm and hence 
ms ~ 0.3 (equation (4)). This, of course, gives us Le and hence Ir = L/Le from 
observational estimates of the bolomctric luminosity, L. We use L/Lx = 10 to 
estimate L from the X-ray luminosity Lx (sect. 4.2) giving 1^ in the range .1 
to .01. In fact, with the radial extent of the BLR fixed, rather than allowed 
to vary according to equation (4), order of magnitude changes in Ir affect the 
line profile of CIV only through small changes in the core. Our chosen range 
of Ir does, howewer, allow a close match to the line ratios (sect. 4.3) and is 
consistent with our estimate of the numerical factor in equation (4). Finally, 
the injection velocity Vi is fixed by the separation between the blue peak and 
the central peak. 

Fig. 2 corresponds to the following choice of parameters: black hole mass 
TOg = 0.3, corresponding to lr{max) =0.1; flux of energy in the nuclear wind 
Lm/LE = 5; normal component of the cloud injection velocity = 2 x lO'^ cm 
s~^. In addition, only clouds with column densities corresponding to 10^'^ cm~^ 
at the inner radius are considered. 

The continuum has been taken as the Mathews and Ferland spectrum for 
this figure, although a power law gives essentially the same result. The cloud 
injection velocity is consistent with a Compton temperature rather higher than 
that for a Mathews and Ferland spectrum, but might be appropriate for a power 
law with no EUV excess (Sect. 1). In any case much of the radiation incident 
on the disc surface will be that emerging from the backs of the clouds forming 
close to the surface and will consequently be in the form of hard X-rays with a 
raised Compton temperature. We therefore believe the required high injection 
velocity to be reasonable. 

In our previous work the fit to the CIV line was made with the inner and 
outer radii of the BLR, the cloud density and the disc flaring chosen freely. 
These arc now, in principle, derived quantities. We have allowed small factors 
between the order of magnitude estimates of these quantities we have made 
for the model in section 2 and the values chosen for the figure. We have used 
Tin = 1.4 X 10^^' cm and Tout — 1-0 X 10^^ cm for the inner and outer radii of 
the BLR, respectively, and n = 3.6 x 10^^ cm~^ for the density at ri„. The disc 
flares by 30° (sect. 2.7). 

The best fit for the line of sight, given the parameter values above, is 58° to 
the disc normal. This fit we have judged by eye: it can be seen from figure 2 
that the smooth model curve cannot fit the small scale detail of the profile, a 
problem which is apparent in the x^-values. In judging our claim that the model 
does fit this observation one should bear in mind the ciifficudticis that all models 
have in reproducing even the broad features of the profiles for this particular 
galaxy. We should not expect the sort of fine detail in figure 2 to emerge here; 
we might expect it to be associated with the fact that clouds are injected into 



14 



the wind in bunches. 

Our line of sight passes close to the disc surface as required if the clouds are 
to partially cover the X-ray source. In addition, the light curves arc most simply 
explained if there is emitting cloud material in the line of sight. The absence 
of a UV bump in this galaxy might be accounted for by sufficiently anisotropic 
emission of the UV, although it is difficult to see how this might be achieved 
physically. 

Note that the double peakiness of the computed line in Fig 2 is a result of 

a complex emission profile, and docs not involve absorption. However the dip 
between the peaks is shallower in the model than in the observed line and almost 
absent if the disc albedo is non-zero (Fig 2b). This suggests that some true 
absorption must be involved in this feature. Such absorption may arise if the 
inner CIV emitting clouds are obscured by clouds further out. This is possible 
in NGC 4151 because of the special angle to the line of sight. The variability in 
the depth and position of the absorption is also in principle compatible with this 
picture. The details do not affect the parameter values we obtain from the profile 
fitting and will be discussed elsewhere. We conclude that the observed CIV line 
profile constrains the model parameters quite tightly and that the predicted 
values of derived quantities provide significant tests of the self-consistency of 
the model. 

The disc itself is assumed to obscure either partially or totally the line emis- 
sion from clouds on the bottom side. (In section 4.6 we shall argue in favour of 
total obscuration and a non-zero albedo, but we need to investigate the prop- 
erties of the partially obscured model to do this.) Without any obscuration of 
clouds below the disc the line would be triple-peaked. In the case of partial 
obscuration the line emission from each cloud is assumed to be blocked wher- 
ever there is material moving in the disc along the line of sight with a relative 
velocity that is within a thermal width of the line centre. What appears as the 
red peak in the observed line is therefore the central peak of the unobscured 
profile. This will be relevant to the behaviour of the two peaks when wc con- 
sider variability in section 4. Finally, in section 4 wc shall be led to consider a 
component refiected from the disc surface. In figure 2 we show the fit for both 
a partially obscuring disc with an albedo of zero and for a totally obscuring 
disc with an albedo of 0.3. (In the latter case the disc is, of course, taken to be 
optically thick out to large radii, so the optically thin and partially obscuring 
regions in Fig. 1 are absent.) 

4 Comparison With Observations 

Detailed studies of NGC 4151 have been prosecuted now for about 25 years. We 
attempt to summarise here what appear to be the main results of these obser- 
vations. Given the determination of the model parameters by the observations 
presented so far, we consider how the remaining data might be interpreted. 
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4.1 Luminosity States 

The profile fitting of CIV fixes the inner radius of the BLR (hence the black 
hole mass) and the properties of the wind. In intermediate luminosity states the 
narrower CIV and Mgll profiles arc obtained by changes in the wind parameters 
associated with decreased luminosity (Sect. 4.2). At low states most of the BLR 
must disappear. This could be attributed to obscuration of the inner clouds or 
to the suppression of cloud formation in the BLR. Obscuration would have to be 
restricted to the clouds (and not affect the central continuum source) since there 
is no correlation with, for example, X-ray absorption features. It might therefore 
refer to an obscuring disc; in any case it would require a special line of sight, 
not only for NGC 4151 but presumably also for all similar intermediate Seyferts 
and for radio galaxies that lose their broad lines intermittently (for example, 
3C390.3). In NGC 5548 there was no evidence for obscuring material in X-rays 
when the broad lines disappeared (Loska et al. 1993). We therefore propose 
to associate the low state with a lack of inner region clouds. A unique feature 
of the model is the short timcscalcs over which the inner clouds are removed 
once injection is stopped (of order a cloud sound crossing time, SxlO^rJg s). 
Thus, the time lag for the change between intermediate Seyfert types following 
a transition in the luminosity state is by the light-crossing time of the inner part 
of the BLR, not by structural changes in the BLR or of some outer obscuring 
regions. 

4.2 Line Widths 

In its relatively infrequent highest luminosity state the relative widths of the 
CIV, Mgll and CIII] lines are exceptional amongst AGN. At their broadest 
the FWHM for CIV and Mg II are between 4000 and 6000 km s^i. The cor- 
responding FWZI are 30 000 and 20 000 km s'^. Unusually for AGN, CIII] 
A1909 is much narrower with a FWHM of 1600 km s"^ and FWZI of 3500 km 
(Clavel et al 1990). This arises from the high density of the inner BLR 
clouds (> 10i°cm-3) which supresses the CIII]A1909/CIVA1549 ratio. Some- 
times Mgll shows the characteristic double peak of CIV and sometimes the blue 
peak is either absent or very weak (Ulrich et al. 1984, 1991). In the low state 
the CIV and Mgll lines are much narrower, around 4000 km s~^ HWZI, with 
suggestions of a broad base, while CIII] is unchanged. Thus, in the low state, 
all three lines have similar profiles. 

The model fit to the profiles requires both a small radial extent for the BLR 
and high density clouds. This results in the similarity of Mgll and CIV line 
widths (FWHM) while CIII] is much narrower. Although this is not shown, Mgll 
is somewhat narrower than CIV, as required, with the difference diminishing at 
lower luminosity. 

However, Mgll does not always show the characteristic double peak of the 
CIV line. Sometimes the blue peak is either absent or very weak (Ulrich et al. 
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1991). We can obtain such behaviour by supressing the Mgll emission from 
clouds immediately leaving the disc which also gives a narrower line profile in 
agreement with Ulrich et al. However, such behaviour does not appear to arise 
naturally in the model. 

In Fig. 3 we show the effect of changing the obscuration properties of the disc 
compared with the CIV line profile observations on Mar 21 and Mar 25 1991 from 
Ulrich et al. (1991). We allow the disc to become completely opaque at less than 
1 X 10^^ cm and at less than 5 x 10^^ cm respectively, so emission from the cloud 
distribution on the far side of the disc is more completely removed. Equivalently, 
we may regard these figures as showing the effect of a changing albedo for a 
completely obscuring disc. We have not been able to mimic these changes by 
any other combinations of alternative parameters or as variability effects (which 
would be included as a radial dependence in the cloud illumination) . The change 
on the red side of the blue peak (with little or no change on the blue side) is 
a signature of a flared disc, the inner clouds initially moving inwards as they 
leave the disc. 

We show the effect on the CIV line of a reduction in the efficiency of cloud 
acceleration in fig. 4. The width of the line is reduced, although this reduction 
is limited by the rotational and injection components of the cloud motions. 
Penston et al. (1981) show a sequence of CIV profiles between Feb 1978 and 
Jan 1979 as the continuum varied from the high to intermediate state. The 
most prominent feature is the disappearance of the profile wings. We would 
therefore like to associate the reduction in cloud acceleration with the decrease 
in luminosity. This cannot be the direct effect of radiation pressure acting on 
the clouds which cannot accelerate clouds to more than the internal sound speed 
over a sound crossing time. Since the survival time of the inner clouds (of order 
lO^s) is less than the lO^s it takes for the nuclear wind to reach the BLR it is 
possible for a reduction in the wind power to influence the cloud velocities in 
the wings of the lines. Since the cloud velocity depends on Lm/ Lv^) (equation 
(psi)) the line width will decrease with L if the flux of momentum in the wind 
decreases with L more strongly than direct proportionality. If Lmlvw oc 
then fig. 3 represents the effect of an order of magnitude change in luminosity 
after the effect of a reduction in wind power has had time to propagate through 
the BLR.. 

Further reduction of the cloud acceleration has no effect on the profile widths 
since the rotational velocity now dominates. The resulting profiles are too wide 
to represent the low state observations. However, as the luminosity drops the 
inner radius of the BLR will move out (equation (4)). We obtain the further 
reduction in the profile width required in the low state if we simply suppress 
the cloud formation within a radius of about 6.5 x 10^^ cm, consistent with our 
estimate for Ir = 0.01 of 5 x lO^^cm. This is in agreement with Peterson (1988) 
and Ulrich (1986) who show that in the lowest state virtually all of the line 
emission arises from outside the high density region in the BLR. The remaining 
BLR clouds provide a hne of the observed FWHM (Fahey et al. 1991). Even 
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so there is still a marked deficit in the core; in our model this missing flux must 
come from a distinct region. Fig. 5 shows the comparison with the data of 
Fahey ct al. (1991). The width of this missing emission is similar to that of the 
CIII] line and the ratio of flux in the CIV core to that in the broad component 
of CIII] is normal for AGN! 

Our profile for CIII] from the BLR, normalised to the same peak intensity, 
would be similar to CIV, in complete contradiction to that observed. However, 
the high inner cloud density means that the flux in this component of the CIII] 
line is merely 10 per cent of the observed flux. The remaining flux consists of 70 
per cent from the NLR (Ulrich et al. 1984) leaving a missing 20 per cent which 
we attribute to the ILR. The computed BLR profile then forms a shallow broad 
base to the observed line. 

These observations clearly suggest that the bulk of the CIII] is coming from a 
separate region in which the CIV/CIII] ratio is typical. This is the intermediate 
line region (called BLR3 in Ulrich et al. 1986). It is interesting to note that 
the ILR is a prediction of the model resulting from just the fit to the profile 
of CIV in the high and low states. We get additional confirmation of this 
from another set of low state profiles (Ulrich et al. 1984) in which the BLR 
contribution is almost entirely absent and the CIV, Mgll and CIII] profiles are 
very similar. We can show that the thermal disc wind gives rise to sufficient 
material to account for the CIII] emission. The flux in CIV is of order 10"'^"'^ erg 
cm^^ s^^ which translates (sec Tabic) to a luminosity in 11(3 of order IO^Lq 
at 2GMpc. (We have not found a directly measured value for the 11(3 flux.) 
Since CIII]/H/3 ^ 3 whereas photoionisation models of the ILR type clouds 
typically give CIII]:H/3 1, about 1/3 of the H/3 flux must be from the ILR. 
From standard recombination theory we can estimate the mass in the ILR as 
30 X M0(L(_ff/3)/lO^L0)(lO^cm~"^/ne). The average density is in the range 
10* cm~^ to 3 X 10^ cm""^ between rc and 3rc (equation (5)). Thus the mass 
required to provide all the CIII] emission is between 1-10 Mq. This is close 
to our estimate of the mass available (equation (30)). Thus the disc wind 
can provide sufficient material to account for the CIII] emission from the ILR, 
although we cannot rule out some contribution from other sources. Furthermore 
the existence of the ILR is reinforced by variability observations: Ulrich et al. 
(1984) found that CIII] does not vary on a time scale of less than a year. Snijders 
(1990) has detected an underlying weak broad variable component which could 
be attributed to the BLR. This implies that the BLR, as proposed here, and 
the ILR cannot merge but that there is a gap between the outer the BLR and 
the inner ILR. 

It is important to note that NGC 4151 is far from unique in possessing an 
ILR. lUE spectra of 3C 390.3 between 1978 and 1986 (Clavel and Wamstcker 
1987) show changes between broad high state BLR profiles and narrow low state 
profiles which we can attribute almost entirely to an ILR. Yee and Oke (1981) 
compared the luminosity of the nuclear component with that in the Balmer 
lines for 3C 390.3 and 3C 382 and found that the permitted line flux departs 
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significantly from a direct proportional dependence on the nuclear continuum. 
For both objects the 1980 scans show a different behaviour from those of earlier 
epochs. Firstly, both objects showed marked changes in the Balmcr line profiles. 
Secondly, between 1979 and 1980 the Balmer line fluxes of each had changed sig- 
nificantly without an accompanying change in the continuum luminosity. (The 
same behaviour is exhibited by 3C 120 (Oke et al. 1980).) A simple time lag 
cannot explain completely the data for 3C 390.3 because it would produce much 
more scatter in the relation between the line flux and nuclear continuum. There- 
fore Yee and Oke suggested a two component BLR, one component reacting on 
a short timescale and the other on a long timescale. This second component is 
equivalent to our ILR. 

Crenshaw, Peterson and Wagner (1988) foimd that the decomposition of 
profiles into broad and narrow components in 3C 445 leaves a residual compo- 
nent for H/3 and Ha. They claim that this component probably arises from a 
region kinematically distinct from the narrow [OUT] emitting region. 

From the [OIII] profiles of Seyfert 1 galaxies van Groningen & de Bryun 
(1989) found that 10/12 objects possess broad wings implying densities in the 
region of 10® cm~^. They call this a transition region intermediate between the 
BLR and NLR. 

4.3 Line ratios 

The average CIV A1549 / CIII] A1909 and CIVA1549/MgIIA2798 line ratios are 
much higher than is typical in AGN or quasars: Clavel et al. (1990) found ratios 
of 7.2 and 6.3 respectively. However, these values are highly variable and were 
found to correlate positively with the continuum flux. In the December 1989 
outburst Clavel et al. found CIV/CIII] = 11.0 and CIV/Mgll = 8.0. In the 
wings of the lines the ratios were even higher with CIV/ CIII] > 50 and CIV/Mgll 
~ 10. Similarly, Uhich et al. (1991) find CIV/CIII] > 10 for \v\ > 2400 km s"! 
and CIV/CIII] > 20 for \v\ > 4000 km s~^, in this case independent of the value 
of /a(1455A). Nevertheless, in the line cores, the CIV/CIII] ratio is ~ 5, similar 
to values found in typical quasar spectra. 

In table I we show the observed line ratios and the model predictions for 
two representative ionising continua. The Mathews and Ferland spectrum is an 
average over observed AGN. However, NGC 4151 is not an average object. In our 
picture this is simply a consequence of a rather special line of sight. But, even if 
the continuum emission from AGN is anisotropic, in NGC 4151 we see just that 
ionising spectrum incident on the BLR clouds (as modified by the presence of the 
clouds in our line of sight). It is possible therefore that the ionising continuum 
more closely resembles the unreprocessed power law with slope -0.9 (e.g. Nandra 
& Pounds 1994). To illustrate the difference this might make we have shown in 
table 1 the results from the model of a power law slope -1 for various luminosity 
states. (The slope in X-rays is only weakly dependent on luminosity (Yaqoob 
& Warwick 1991).) Also shown for comparison is the photoionisation model 
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calculation of Ferland and Mushotsky (1982) for a homogeneous BLR and an 
extrapolated spectrum corrected for X-ray absorption. It is interesting to note 
the improved agreement for Mgll with the power law spectrum (although this 
conclusion is sensitive to any change in the column density). 

In the high state the observed low GUI] /CIV ratio arises from the high 
density clouds in the BLR. At lower luminosities the contribution from the ILR 
is increasingly important and the ratios are predicted to become closer to typical 
values. The low state ratios do not appear to be given in the literature, but 
since the low state profiles involve the line cores where the high state ratios are 
more normal, we expect this prediction to be borne out. 

Note that, as for other inhomogeneous models, there is no incompatibility 
between the ionization parameter deduced from line ratios and that deduced 
from reverberation maps of the size of the BLR. Note also that for the parame- 
ters of the model appropriate to NGC 4151 the clouds are found to move only a 
relatively small radial distance from their point of origin on the disc; the mean 
gas density in the clouds therefore follows fairly closely n(r) ex 1/r^. (This is 
not the case, in general, in systems where the normal injection speed is lower 
or where the BLR has a larger radial extent.) 

4.4 Blue shifted 'absorption' 

A double peak in the CIVA1549 line was first noted by Anderson and Kraft 
(1969). Cromwell and Wcyman (1970) and Ferland and Mushotzsky (1982) also 
report a variable feature on the blue side of the CIV line which they discuss 
in terms of variable absorption. Profiles showing this feature in a number of 
different luminosity states of NGC 4151 arc presented by Pension et al. (1981), 
Ulrich et al. (1984) and Stoner & Ptak (1986). From these it is clear that the 
absorption feature varies, both in wavelength and relative depth in response to 
the variations in the continuum 

In the model a double peak in the CIV line profile is the result of a con- 
tribution from rotational motion together with a blue shifted component from 
the injection velocity. The corresponding red shifted peak from clouds on the 
far side of the disc is obscured. The position of the trough in the profile moves 
towards lower velocity as the continuum luminosity decreases in agreement with 
observation. Nevertheless the depth of the trough in the model profile in the 
high luminosity state is significantly shallower than observed. It is likely there- 
fore that some true absorption is also occuring in this velocity range, perhaps 
as a result of obscuration of the inner broad line clouds by those at larger radii. 
The overall absorption line structure is complex and we plan to return to the 
subject elsewhere. 
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4.5 Profile asymmetries 

In 1986-88 all the broad lines were asymmetric with more flux in the red wings. 
The mid-points of the FWHMs were shifted blueward with respect to the narrow 
line centres by 1000 km s~^. The narrow lines themselves show a small blue 
shift with respect to the galaxy. Ulrich et al. (1984) find for CIV that the HWZI 
in the blue exceeds that in the red. Yet the high luminosity profile published 
by Fahcy et al. (1991) clearly shows excess emission and larger width in the red 
wing. For Mg II the red and blue wings have similar widths with more emission 
sometimes in the red. Note that the red wing in CIV may be blended with Hell 
A1640 and OIII] A1663 (Clavcl et al. 1987); the Mgll blue wing may be affected 
by He I A2733 or the Fe II multiplets UV62 and 63. Large changes in the ratio 
of intensities in the wings between campaigns and different responses of the red 
and blue wings to continuum variations have been observed by Ulrich ct al. 
(1991). With no reflection or obscuration from the disc the red and blue wings 
of the high state model profile are symmetric, but become red asymmetric with 
the inclusion of a non-zero albedo. Changes in the asymmetry arc obtained in 
the model from changes in the disc albedo, since this affects mainly the intensity 
and lag of the red component. In the low state the model CIV (or Mgll) profile 
is asymmetric towards the blue i.e. the blue wing is wider and has more flux 
(fig. 5). We need to resolve the problem of blending and to include non-steady 
state profiles in order to make a firm comparison with the observations. 

Ulrich et al (1991) also note a transient red component in CIV possibly 
connected with the narrow so-called Li and L2 components which are visible 
only in the lower luminosity states. The Li and L2 features cannot come from 
our BLR. They may arise from gas associated with the linear radio structure. 
The intensities of Li and L2 are not correlated with the continuum flux, but are 
well correlated one with the other. The features have a large velocity separation 
but any time delay between the variations has to be less than 4 days (Clavel et 
al. 1987). Thus a structure close to the plane of the sky is indicated. If this is 
the radio structure then it is consistent with our model. 

4.6 VariabiHty 

Ulrich et al (1984) found typical 2-folding timescales for CIV of 19-47 days and 
for Mg II of 27-72 days; the shortest two-folding timescale over a period of about 
2 years of observation was 7 days for CIV and 20 days for Mgll. The light curves 
of CIV and Mg II are similar to that of the continuum (Clavel et al. 1990) with 
CIV undergoing changes greater than 60 per cent but Mgll less than 50 per 
cent. However, the CIII] line remains constant. The transfer functions for Mgll 
and CIV are similar and show zero delay. This is contrary to Clavel et al. (1990) 
who observed a delay of 4 ±3 days, but any discrepancies can be explained by 
varying continuum pulse lengths. 

The peak of the line response to a continuum pulse gives a characteristic 
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lag time. For NGC 4151 there is effectively zero delay in the initial response of 
the line as a result of the presence of material in the line of sight. For a short 
continuum pulse the peak occurs near zero delay because there is a significant 
fraction of the clouds within a sampling interval of the line of sight. In its low 
state, where the inner BLR clouds are absent, clouds in the line of sight still 
give an immediate response, but the peak response now occurs between 10-20 
days, corresponding to the time at which the equal delay surfaces contain the 
maximum number of clouds. 

For a non-sphcrical BLR, such as we have in the model, this lag depends on 
the length of the continuum pulse as well as the BLR geometry. For a short 
continuum pulse, duration less than or equal to the sampling time, the inner 
clouds dominate the transfer function so the peak occurs at small time delay. 
For longer pulses the delay is determined by the peak of the emission reponse 
at larger radii. Such a dependence has been confirmed for NGC 5548 by Netzer 
(1991). 

In Fig. 6 we show the model response of the CIV line to a delta continuum 
pulse from zero background for a partially obscuring disc with zero albedo. The 
lag is ^ 1/2 day. Mg II and the contribution of the BLR to H/? are predicted to 
respond similarly. The sampling time in the only available observations of the 
transfer function of NGC 4151 (Netzer 1991) is ^ 3 days. Within this resolution 
the line response of CIV is indeed immediate. In the model response there is 
a second peak which is the signature of the Keplerian motion of the clouds; 
there is no evidence for this in the observed data, but the resolution is almost 
certainly too low. 

If we consider the transfer functions for the line wings and core separately 
we find that changes to the bulk of the core lag those in the wings and that, 
for the partially obscuring disc, the blue wing leads the red. The former agrees 
with observation; the latter does not. The 2D transfer function is difficult to 
construct observationally and none appears to be available for NGC 4151. Fig. 
7 shows our prediction for the 2D transfer function for CIV with the same 
conditions as in Fig. 6 and with zero albedo for the disc. It is again clear that 
changes in the blue wing lead those in the red in probable disagreement with 
observation. 

We remedy this discrepancy by taking a non-zero disc albedo. The transfer 
function for this reflected component is shown in fig. 8 for a reflection efficiency 
of 0.3. It can be seen that incorporating reflection will give a transfer function 
in which the red and blue components of the line respond with the same lag to 
continuum changes. 

4.7 Profile variations 

Gaskell (1988) found that variations in the red wings lead the blue in CIV and 
Mg II by a few days. The precedence of the red wing is not confirmed by other 
observers. As Gaskell points out, the exclusion of the low state of April 21 1980 
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in the cross-correlation function changes his result. Clavel et al. (1990) found 
equal delays in the red and blue wings of 3.2±3 days. This is consistent with 
Ulrich et al. (1991) who find delays in the response of the wings less than their 
sampling time of 5 days, although the red and the blue wings do not always 
respond to continuum changes in the same way. Fahey et al. (1991) found that 
the blue peak of the CIV line exhibits stronger variability than the red peak. 
The emission of the red peak is also strongly correlated with the total CIV flux, 
whereas the blue peak seems to be more sensitive to variations in the CIV flux 
and exhibits a sharp turn-on near fciv = 10""'^'^ erg cm~^ s^^. 

Clavel et al. (1990) observed a fractional variation in the flux in the CIV 
wings ( \v\ > 3000 km s~^) of 0.34, much larger than that of the line core 
(|f I <3000 km s"'^), of 0.12. The same pattern is observed in the Mgll lines, 
although it is less pronounced than in the CIV line. Clavel (1991) found that 
the fractional variation of the wing flux is as large as that of the continuum, 
whereas the flux in the core only changes by ^ 50%. 

Changes in the model profile which affect only the red side of the line or the 
red side before the blue can only be a result of changes in the transmittance 
or the reflectance of the disc. In general, if there are no structural changes we 
expect continuum variations to produce simultaneous changes in the red and 
blue wings, if we include reflection from the disc, and for changes on the blue 
side to precede those on the red if there is no reflection. Since the unreflcctcd 
line is asymmetric and the reflected component contributes a higher fraction to 
the red than the blue the relative changes in the wings can be complex. The 
fractional changes in the wings will also depend on the continuum pulse length, 
because of the extra path length for the reflected component from the outer 
clouds. 

Fig. 6 shows the transfer function for the CIV line as a whole. We have 
also constructed a diagram for the core of the line (±3000 km s^^) only. The 
variation in the core provides between .03 and .1 of the peak throughput; the 
remainder is provided by the wings so the fractional change in the wings must 
exceed that in the core. 

The core of CIV responds to UV variations with a delay of 6.5±3 days, 
significantly longer than the response time of the wings which is 3.2±3 days on 
average (Clavel 1991). 

The Id model transfer functions for the wings and line core separately show 
that the bulk of the core lags the wings in response to UV continuum changes. 

4.8 Optical lines 

Antonucci & Cohen (1983) observed that part of the H/3 line flux responded 
to continuum changes immediately, but a large part responded only gradually. 
Most or all of the flux of the higher-series Balmer lines and the A4686 line of 
Hell responded immediately. They observed no gross profile changes. Also, 
the features of the H/3 profile are very similar in 1981 and 1975, through high 
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and low states. In the first half of 1984, while the nucleus was in a low state, 
De Robertis (1985) observed little, if any, broad H/3 and very weak broad Ha 
emission. 

Comparisons between H/3 and CIV show that CIV is much broader than 
H/3 (Osterbrock & Koski 1976); this is especially evident in the low luminosity 
state. The change in H/3 between the high state and the low state was found to 
be qualitatively similar to that of the CIV profile (Ptak & Stoner 1985). The 
Lya/H/3 ratio is of order 11.6 (Ferland & Mushotzky, 1982), much higher than 
typical for active galaxies; the intrinsic value could be even higher since the 
blue side of the Lya line is missing (Penston et al. 1979, Ferland & Mushotzky 
1982). Boksenberg & Shortridge (1975) observed that the Hel AA5876 and 7065 
lines arc twice the width of the Balmer line wings. 

Observations by Lyutyi (1976) show that the Ha + [Nil] emission feature 
varies in intensity with a 30-d lag for the 1970-72 period and 60-70 days in 
1973-76. Maoz et al. (1991) have observed a lag in the Balmer lines of 9±2 
days, still significantly longer than the UV line lags. 

For the dense inner BLR clouds with our high ionisation parameter we find 
that the Balmer lines, are suppressed relative to CIV (For Ha this is clear from 
Table 1.) Thus the ILR is a significant contributor to the Balmer lines. This is 
consistent with the relative line widths, ratios and time lags. At lower luminosity 
the BLR contribution is relatively less important for the CIV line, reducing the 
CIV/H/3 ratio and giving lines that differ less in width. 

4.9 Continuum 

Clavel et al. (1990) found that in the UV the typical amplitude of variability 
was the same at 1715A and 1455A but there did appear to be some spectral 
variations which preceded the changes in flux levels by about 12 days. The ratio 
of maximum to minimum flux was 3.1 for the continuum at 1455A but only 1.24 
at 5000A. The continuum below 1900A underwent small but significant spectral 
variations. (Clavel et al. 1990). The X-ray luminosity in the 2-10 keV band has 
been found to vary between ~ 2 — 20 x 10''^ erg s^^ (Yaqoob et al. 1993). The 
spectral variations are complex but consistent with a picture of an intrinsically 
variable source, with a spectrum that steepens on brightening. This is seen 
through an absorbing column that has been modelled as cold (T < lO^K) gas 
that partially covers the source and varies on time scales of months to years 
within a range of column densities Nh ~ 5 — 15 x 10^^ cm~^ (Barr et al. 1977, 
Pounds et al. 1986, Yaqoob et al. 1989) or as warm gas where the soft X-ray 
opacity is reduced by partial ionization (Warwick & Done 1994, see also Weaver 
et al. 1994). The depth of the iron K-edge implies an iron overabundance of 
order twice solar. There is very little reddening of the optical continuum with 
A„ « 0.2 (Mushotzky et al. 1978). Therefore the X-ray absorption is not due 
to dust. 

We have taken the spectral shape to be independent of luminosity state. In 
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fact, modest changes in spectrum do not appear to alter qualitatively any of the 
conclusions of the paper. 

We note the possible connection between our BLR model and the X-ray 
observations. From fig. 1 we see that if the X-ray source is located centrally 
then there will be some absorption from the ILR and BLR clouds in the line of 
sight. The ILR is a candidate for the fixed component of the absorbing screen 
(Yaqoob et al. 1989) with the BLR clouds forming a variable partial coverer. 
Yaqoob et al. (1991) excluded the warm absorber model for this source, although 
Weaver et al. (1994) include a warm absorber in a two component picture. In 
our model the X-ray source must be partially covered but there could be also a 
scattered component and possibly some absorption from the warm cloud debris. 
The lack of reflection features from a disc (Maisack and Yaqoob 1991) is another 
indication that the angle to the line of sight is closer to edge-on than face-on. 
The complex absorption cut-off observed by Ginga is also difficult to understand 
unless the line of sight is closer to being edge-on than face-on (Yaqoob 1993). 
This provides some support for the geometry of our model. 

5 Conclusions 

We began by constructing a theoretical model for the injection and dynamics 
of clouds in the broad and intermediate emission line regions. In this model 
the BLR clouds arise from the interaction between a nuclear wind and the 
radiatively heated disc surface and the ILR clouds form from the interaction of 
the wind with a thermally driven outflow from the disc. When applied to NGC 
4151 we found we could specify a complete set of parameters fairly tightly from 
the profile of the CIV A 1549 line in the high luminosity state of the nucleus. In 
particular the density of the inner clouds is high (~ 10^^ cm~^) and the region 
is small (10^^ cm < r < 10^^ cm). From this and the assumption that the power 
in the wind drops strongly with decreasing radiative power, we derived profile 
fits to CIV in some intermediate luminosity states and accounted for the profile 
widths of other UV lines and the Balmer lines in various states. From the low 
luminosity state behaviour we deduced the existence of an extended intermediate 
line region with cloud densities < 10^ cm~^. This picture agrees with that 
derived from the transfer function and correlation analysis of the response of 
the lines to continuum variability. Provided we include reflection from the disc 
in constructing the line response then the model is also in qualitative agreement 
with what is known of the transfer function (the relative responses of red and 
blue wings and line core). The geometry deduced from our model BLR also 
leads us to suggest an interpretation of the X-ray absorption data to which we 
shall return in more detail elsewhere. 

Of other possible models for NGC 4151 pure disc motion appears to be ruled 
out since variations in the line core would then precede that in the wings. In the 
Mathews model of oscillating clouds (Mathews 1993) the reponse of the core 
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and wings should be simultaneous, which it is probably not. In the picture of 
Dyson and Perry (1985) it is difficult to see how the red and blue wings of the 
line profiles could change simultaneously in a spherical outflow. A biconical 
model could be made to fit the available data for NGC 4151, but this would 
probably require us to treat this galaxy as an exception since the 2D transfer 
function of NGC 5548 (Horne et al. 1991) appears to exclude a biconical model 
for that galaxy. 

In our model NGC 4151 is a 'special' case only in as much as we require a 

particular orientation to the line of sight. Wc shall show elsewhere that essen- 
tially the same model can be constructed for NGC 5548, except that the line of 
sight is at greater inclination to the disc. We shall also show how NGC 4151, 
which is usually regarded as an exceptional case, can be fitted into a unified 
picture for active galaxies. 
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Figure Captions 



Fig. 1 Sketch of the line emission regions of NGC 4151. The BLR hes between 
and rout in a region of the disc dominated by external illumination where 
strong Kelvin-Helmholtz instabilities at the interface with the nuclear wind, 
labelled v^, lead to cloud formation and entrainment. The ILR lies between 
ric where Compton heating is sufficient to drive material from the disc and 
riso where the heating timescale exceeds the escape time. The disc is probably 
optically thick at small radii and thin at large radii; various assumptions are 
explored in the text. 

Fig. 2. The observed high state CIV profile (solid line) compared with the 
model profile for a disc of zero albedo which partially covers clouds on the far 
side (dotted line) and for a disc with an albedo of 0.3 with clouds on the far 
side totally obscured (dashed line). The model lines have been scaled vertically 
to give the best fit. Velocities on the horizontal axis are given in km s~^ with 
negative velocities to the red. 

Fig. 3. The effect of changing disc obscuration. The solid line is for a disc 
which obscures totally at distances < 5 x 10^® cm and the dashed line, on the 
same vertical scale, is for total obscuration at < 1 x 10^^ cm. The dotted line is 
the profile difference. From the symmetry of the model these profiles effectively 
show also the effect of increasing albedo (from to 1 in the region 1 x lO^^cm 
to 5 X lO^'^cm). 

Fig. 4. A sequence of 3 CIV model profiles, on the same vertical scale, from 
high to intermediate states. The profile narrows as the acceleration of clouds 
by the wind decreases corresponding to decreasing luminosity states. 

Fig. 5 CIV profile fit for the low luminosity state. The observed profile is 
shown as the solid line, the BLR contribution is the dot-dashed line, scaled 

vertically to fit the bulk of the line outside the core, and the difference is the 
dashed line. The difference profile is the contribution from the ILR in the model. 

Fig. 6 One dimensional model transfer function for the CIV line (a) in a 
low-luminosity state {Ir = 0.01, dashed line) and (b) in a high state {Ir = 0.1, 
solid line). The number of half day intervals is marked on the x-axis. 

Fig. 7 Two dimensional model transfer function for the CIV line (with a 
partially obscuring disc and zero albedo). The time is given in half-day inter- 
vals and the velocities in km s~^. 

Fig. 8 Two dimensional model transfer function for the component of the 
CIV line reflected by a disc with constant albedo. The time is given in half-day 
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intervals and the velocities in km 



